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Abstract 

We study the stability of the two-neutrino vacuum oscillation solution of the solar neutrino 
problem with respect to changes of the total fluxes of 8 B and Be neutrinos, $b an d ^Bc- 
For any value of <E>Bc from the interval 0.7$g^ < <3?bc < 1.3<&§ P the solar v e oscillations 
into an active neutrino, v e «-> ^ M ( T ), provide at 95% C.L. a description of the existing solar 
neutrino data for <1>b — (0.35 — 3.4) ^jp, and being the fluxes in the solar model 
of Bahcall-Pinsonneault from 1992. For 4» Bc = (0.7 - 1.3)4>| P we find also at 95% C.L. two 
new (one new) v e <-> *V T ) [y e <-> u s ) oscillation solutions: i) for $b — (0.35 — 0.43)4>g P at 
Am 2 ^ (4.7-6.5) x 10~ 12 eV 2 ((4.8-6.4) x lO" 12 eV 2 ) and sin 2 26 k, 0.71 (0.74), and ii) for 
$B = (0.45 - 0.65)$| p at Am 2 ^ (3.2 - 4.0) x 10" 11 eV 2 and sin 2 26 ^ 0.59. The physical 
implications of the new solutions for the future solar neutrino experiments are discussed. 
The data rule out at 97% - 98% (99 %) C.L. the possibility of a universal (neutrino energy 
independent) suppression of the different components of the solar neutrino flux, resulting 
from solar v e oscillations or transitions into active (sterile) neutrino. 



*Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 
BG-1784 Sofia, Bulgaria. 



1 



1. INTRODUCTION 

In the present paper we investigate the stability of the vacuum oscillation [1-4] solution 
of the solar neutrino problem [5,6] with respect to variations of the total fluxes of the solar 
8 B and 7 Be neutrinos. Recent studies have indicated that the current solar model predictions 
[7-12] for the 8 B neutrino flux, $ B , vary from model to model with rather large uncertainties 
[12,13]. The results for $ B derived in all solar models presently discussed in the literature 
except that of ref. [12] , lie in the interval (4.43—6.62) x 10 6 v e j cm 2 /sec, while the prediction of 
the "low" flux model of ref. [12], $ B = 2.77x 10 6 z/ e /cm 2 /sec, differs from those of the "high" 
flux models of refs. [7] and [11] approximately by the factors 2.0 and 2.4. The predictions [7- 
12] for the total flux of 7 Be neutrinos, <3> Be , vary by ~20%, from <3> Be = 4.34 x 10 9 z/ e /cm 2 /sec 
in ref. [8] to $ Be = 5.18 x 10 9 z/ e /cm 2 /sec in refs. [11]. At the same time none of the solar 
models developed to date provides a satisfactory description of the existing solar neutrino 
data [5,14-16]. In particular, the upper limits on the value of the 7 Be neutrino flux, which can 
be inferred from the data, are considerably lower than the values predicted by the models, 
as first noticed in ref. [17] and confirmed in several subsequent more detailed studies [18] 
utilizing different methods. The above result follows not only from joint analyses of the data 
from all solar neutrino experiments, Homestake [5], Kamiokande [14], GALLEX [15] and 
SAGE [16], but also from the Homestake and Kamiokande, or from the Kamiokande and 
SAGE and/or GALLEX data. Since the recent calibration of the GALLEX detector [19] 
leaves little room for doubts about the correctness of the GALLEX results, both the data 
from the Davis et al. and Kamiokande experiments have to be incorrect in order for the 
indicated conclusion to be not valid. The discrepancy between the value of $ Be suggested 
by the analyses of the available solar neutrino data and the solar model predictions for 
$ Bc represents a major new aspect of the solar neutrino problem. No astrophysical and/or 
nuclear physics explanation of this discrepancy has been proposed so far. 

Assuming that the 7 Be neutrino flux has a value in the interval 0.7$b^ < $ Bc < 1.3$b^, 
where is the flux predicted in the reference solar model of Bahcall - Pinsonneault [7], 
we determine in the present study the range of values of the 8 B neutrino flux, for which the 

2 



results of the solar neutrino experiments can be described in terms of two-neutrino vacuum 
oscillations of the solar neutrinos into an active u e <-> iV T ), or sterile u e <-> z/ s , neutrino. 
Similar analyses for the MS W solution [20] with solar v e transitions into an active neutrino, 
v e — > v ix{t)i were performed in refs. [21,22]. Results for the case of solar v e <-> v^t) oscillations 
were obtained in ref. [23] for $ Bc = 0.8$^ and 0.4$| p < $ B < 1.6<^>| p , where is the 
8 B neutrino flux predicted in the reference model [7]. However, we find, in particular, that 
at 95% C.L. a new vacuum u e <-> v^t) { v e u s) oscillation solution of the solar neutrino 
problem exists in the region Am 2 = (4.7 - 6.5) x 10~ 12 eV 2 ((4.8 - 6.4) x 10~ 12 eV 2 ) and 
0.71 (0.74) £ sin 2 29 < 1.0 for $ B = (0.35 - 0.43)$i P and $ Bc = (0.7 - 1.3)$| P , Am 2 
and sin 2 29 being the two parameters characterizing the oscillations (see, e.g., refs. [1- 
4]). A second v e <-> v^r) oscillation solution is found for 0.45$ B P ^ $ B ^ 0.65$ B P and 
0.7$| P < $ Bc < 1.3$ B P , and for values of Am 2 and sin 2 29 which lie within the intervals 
Am 2 = (3.2-4.0) x 10~ n eV 2 and 0.59 £ sin 2 29 < 1.0. Both these solutions were not noticed 
in ref. [23]. For Am 2 > 4.1 x 10~ n eV 2 the v e <-> oscillations allow to describe at 95% 
C.L. the existing solar neutrino data for any value of $ Be from the interval (0.7 — 1.3)$ PP (for 
$ Bc = 0.7$ B P ) and for 0.57 (0.51)$ B P ^ $ B ^ 3.4$^ p . The corresponding allowed regions 
of values of the parameters Am 2 and sin 2 29 in all these cases are derived as well. Except 
for 0.35$ B P ^ $ B ^ 0.44$ pp the oscillations into a sterile neutrino v e <->■ v s are excluded at 
95% C.L. if $ Bc = (0.7 - 1.3)$g P ; at 98% C.L. they are allowed by the (mean event rate) 
solar neutrino data for 0.32$| p & $ B £ 1.8$| p . 

We have also performed a study which shows that the data do not favour the hypoth- 
esis of neutrino energy independent suppression of the solar neutrino flux: it is excluded, 
depending on the value of $ Bc from the interval (0.7 - 1.3)$| P , at (97% - 98%) C.L. when 
the suppression is due to u e — > v^r) {y e <-> ^( T )) or u e — > transitions (oscillations), and 
at (99.0% - 99.7%) C.L. if it results from u e — > v s {y e <-> v s ) transitions (oscillations). 

The vacuum oscillation solution at Am 2 ^ 4.4 x 10 -11 eV 2 imply a non-negligible sup- 
pression of the pp v e flux (approximately by a factor (0.50 - 0.70)), and a not very strong 
suppression of the 0.862 MeV 7 Be u e flux, the relevant suppression factor ranging from ap- 
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proximately 0.30 (possible for $ Bc = 1.3$|£) to 0.98 (possible if $ Be = (0.7 - 1.0)$j£ and 
$b = $i P ). 

The new z/ e <-> z/ M(r ) {y e <-> i/ s ) oscillation solution located in the region Am 2 = (4.7 - 
6.5) x 10~ 12 eV 2 , which exists for low values of $ B , is quite similar to the MSW nonadiabatic 
[y e — > v^t) [21] or u e — > z/ s transition [26]) solutions for similar values of <3>b: it corresponds 
to a rather strong suppression of the 0.862 MeV 7 Be v e flux (by a factor (0.06 - 0.26)), to a 
moderate suppression of the pp neutrino flux (by a factor not less than 0.7 only) and to 8 B 
neutrino flux practically not affected by the oscillations. The physical implications of the 
indicated new vacuum oscillation solutions for the future solar neutrino experiments are also 
briefly discussed. Our results show, in particular, that the v e <-> v^t) vacuum oscillation 
solution of the solar neutrino problem is stable with respect to changes of the predictions 
for the 8 B and 7 Be neutrino fluxes. 

We use the latest published data from all four solar neutrino experiments [5,14-16] in 
our analysis: 

R(Ar) = (2.55 ± 0.25) SNU, (1) 

$ B xp = (2.89 ± 0.42) x 10 6 cm- 2 sec~ 1 , (2) 

R G ALLEx(Ge) = (77.1 ± 9.9) SNU, (3) 

RsAGE(Ge) = (69 ± 13) SNU, (4) 

where R(Ar), and RcALLEx(Ge) and RgAGE(Ge), are respectively the average rates of 37 Ar 
and 71 Ge production by solar neutrinos observed in the experiments of Davis et al. [5], 
and GALLEX [15] and SAGE [16], and $ B xp is the flux of 8 B neutrinos measured by the 
Kamiokande collaborations [14]. In eqs. (1) - (4) the quoted errors represent the added in 
quadratures statistical (1 s.d.) and systematical errors. 

2. THE 8 B AND 7 Be NEUTRINO FLUXES 

It is convenient to introduce the parameters 



in terms of which we shall describe the possible deviations of $ B and $ Be from their values 
in the reference model [7]. The fluxes $ B and <£> Be i n the models [7,8,11,12] correspond, 
respectively, to f B = 1.0; 0.78; 1.16; 0.49, and f Be = 1.0; 0.89; 1.06; 0.89. 

The Kamiokande data, evidently, imposes limits on the values $ B (and f B ) can possibly 
have. The expression for the predicted event rate in the Kamiokande detector, R(K), if 
the 8 B (electron) neutrinos undergo two-neutrino transitions into an active neutrino v^m 
(due to vacuum oscillations v e «-> or MSW transitions v e — > i^( T )), or v^r) (due to 
spin-flavour conversion u e — > u^m) on their way to the Earth, has the form Q 

14.4 MeV 

R(K) = f B $i P J n(E) a K (E) [P(E) + 0.16(1 - P(E))] dE, (6) 

14.4 MeV 

where n(E) is the normalized to 1 spectrum of 8 B neutrinos, / n(E)cffi = 1, ctk(E) is the 
v e — e~ elastic scattering cross-section for 8 B neutrinos with energy E, in which the recoil e~ 
detection efficiency and energy resolution functions of the Kamiokande detector are included, 
P(E) is the probability of survival of the 8 B v e having energy E ((1 - P(E)) is the probability 
of the v e — > Vn{f) transition due to vacuum oscillations or the MSW effect, or of the u e — ► v^t) 
conversion), and we have used the fact that a v , )e (E) / a Vee (E) = ap,. e (E)/a l/ee (E) = 0.16 
in the energy range of interest, a Uie (E) and <Tp ie (E), l=e,/i, r, being the v% — e~ and v\ — e~ 
elastic scattering cross-sections. In the case of v e <-» v s oscillations or v e — > v s transitions 
the term with the coefficient 0.16 is absent from the expression in the right hand side of eq. 
(6). 

Given R(K), $g P , n(E) and <tk(E), the minimal allowed value of f B , as it follows from 
(6), is determined by the maximal possible value of [P(E) + 0.16 (1 - P(E))], which is 1 
and is reached when P(E) = 1. Thus, we have f B > R(K)/R B p(K) = (0.51 ± 0.07), where 
R B p(K) is the event rate predicted in the BP model [7], and we have used the Kamiokande 



1 In the numerical calculations we have performed we have included the Kamiokande energy 
resolution and trigger efficiency functions in the expression under the integral in eq. (6). 
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result, eq. (2). At 99.73% (95%) C.L. this implies 



f B £ 0.30 (0.37). 



(7) 



It is trivial to convince oneself that the above lower limit on f B holds also in the case 
of solar v e two-neutrino oscillations (transitions) into sterile neutrino v s , as well as for 
oscillations (transitions) involving more than two neutrinos (sterile and/or active). The 
limit (7) is universal: it does not depend on the type of possible oscillations (transitions), 
and on the specific mechanism responsible for them. 

Similarly, the maximal allowed value of f B by the Kamiokande data corresponds to min 



[P(E) + 0.16 (1 - P(E))] = 0.16. We have then: f B < R(K)/(0.16 R BP (K)) = (3.2 ± 0.44), 
which gives at 99.73% (95%) C.L. 



Inequality (8) is universal for two-neutrino solar v e oscillations or transitions into an active 
neutrino z/ m(t) or v^ r) . 

Contrary to the lower limit (7), the upper limit (8) is not valid for two-neutrino u e <-> v s 
{y e — > u s ) oscillations (transitions) or u e oscillations (transitions) involving more than two 
neutrinos. In the first case, for instance, the maximal value of f B would correspond to the 
min P(E), and the use of the general property of the probability P(E), min P(E) = 0, does 
not allow one to derive a useful upper limit on f B from the Kamiokande data. 

In our study of the stability of the results on the vacuum oscillation solutions with respect 
to $ B and <3> Be variations the following approach is adopted. The fluxes of the pp, pep and 
the CNO neutrinos (see, e.g., refs. [6]) are kept fixed and their values were taken from ref. 
[7]. The fluxes of the 8 B and 7 Be neutrinos, and correspondingly, f B and f Be , are treated as 
fixed parameters, which, however, are allowed to take any values within certain intervals. In 
the case of <3> Be the interval chosen corresponds to 



f B £ 4.5 (4.1). 



(8) 



0.7 < f Be < 1.3. 



(9) 
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It is somewhat wider than the interval formed by the current solar model predictions: 0.89 
- 1.06. For <3> B values in the intervals determined by the inequalities (7) and (8) were 
considered. The searches for a u e <-> u s oscillation solution were preformed for 0.3 < f B < 4.0. 

The indicated approach was motivated by the fact that the contributions of the CNO 
neutrinos to the signals in all three types of detectors [5,14-16] are predicted to be relatively 
small [6-12], and that (apart from the CNO neutrinos) the spreads in the predictions for 
the fluxes $ B and <3>Be are the largest. Some of the values of $Be used in the analyses, as 
those corresponding to fee = 0.7 and 1.3, for example, are incompatible with the constraint 
on the solar neutrino fluxes which the data on the solar luminosity impose (see, e.g., refs. 
[27,28]): 

$ pp + 0.958$ Bc + 0.955$ CNO + 0.910$ pcp = (6.517 ± 0.02) x 10 10 cm^sec -1 , (10) 

where $cno — ^n + ^O; and $ pp , $ pep , $n and $ are the fluxes of the pp, pep and the CNO 
neutrinos. However, a 20% - 30% change in $ Bc with respect to $ B g = 4.89 x 10 9 z/ e /cm 2 /sec 
is required by (10) to be balanced by only a few percent change of the pp neutrino flux, and 
the latter will have a small effect on the predictions for the signal in the Ga-Ge experiments 
[15,16]. Besides, the aim of the present study (as like of the analogous studies of the MSW 
solutions in refs. [21,22,26]) is, in particular, to determine the ranges of values of $ B and 
<3>Be for which the possibility of vacuum oscillations of solar neutrinos cannot be excluded 
by the existing solar neutrino data. Certainly, values of <3> B corresponding to, e.g., f B = 3 
seem at present unlikely to appear in any realistic solar model. 

In the absence of "unconventional behaviour" (vacuum oscillations, MSW transitions, 
etc.) of solar neutrinos, the signals in the Cl-Ar and Ga-Ge experiments can be written in 
the following form within the above approach: 

R(Ar) = (6.20f B + 1.17f Bc + 0.40 CNO + 0.23 pep ) SNU, (11) 

R(Ge) = (70.8 PP + 3.1 pcp + 35.8f Bc + 13.8f B + 7.9 C no) SNU, (12) 
where 6.20f B SNU is the contribution in R(Ar) due to the 8 B neutrinos, etc. 
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We have used the x 2 — method in our analysis. In computing the y 2 for a given pair of 
values of the parameters Am 2 and sin 2 29 we have ignored the estimated uncertainties in the 
reference model predictions [7] for the solar neutrino fluxes as the ranges within which we 
have varied $ B and $ Be exceed by far the uncertainties. We did, however, take into account 
the uncertainties in the detection cross-sections for the detectors [5,14-16]. 

3. THE VACUUM OSCILLATION SOLUTIONS 

3.1 The Case of u e <-> z/ M(r ) Oscillations 

3.1.1 Allowed Regions of the Parameters 

Searching for vacuum u e <-> z/ m(t ) and u e <-> v s oscillation solutions we have scanned the 
region 1(T 12 eV 2 < Am 2 < 10" 9 eV 2 and 10~ 2 < sin 2 29 < 1.0. It was found that at 95% 
C.L. and for 0.7 < fee < 1.3 the two-neutrino v e <-> v^( T ) oscillations of the solar v e allow 
one to describe the data (1) - (4) for rather large intervals of values of f B . These intervals 
depend somewhat on the value of f Be . Below we give the solution intervals for f B (at 95% 
C.L.) in the three representative cases of f Be = 0.7; 1.0; 1.3: 

v e <-> Vn{ T ) '■ fee = 0.7, 0.35 ^ fe ^ 3.4, (13a) 

fee = 1.0, 0.35 £ fe £ 3.4 (136) 

fee = 1.3, 0.35 £ fe £ 0.43 and 0.46 £ f B £ 3.4. (13c) 

The allowed regions of values of Am 2 and sin 2 29 corresponding to the solutions (13a), (13b) 
and (13c) are shown in Figs, la, lb and lc, respectively. 

As Figs, la-lc illustrate, a deviation of the 8 B and 7 Be neutrino fluxes from the values 
corresponding to fe = fe c = 1 leads to two effects: i) noticeable shift (and change in size) 
of the allowed Am 2 — sin 2 29 regions in the case fe = fe c = 1 towards smaller (fe < 1) or 
larger (fe > 1) values of sin 2 29 with the allowed values of Am 2 remaining practically within 
the interval of the fe = f Be = 1 solution, 4.4 x 10 -11 eV 2 £ Am 2 £ 9.8 x 10~ n eV 2 , and ii) 
appearance of new allowed regions of values of Am 2 , i.e. of new solutions, at Am 2 ^ 4.1 x 
10" 11 eV 2 . We find two such new solutions (see Figs, la - lc): 
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(A) for 0.7 < f Bc < 1.3 and 0.35 & f B & 0.44 with Am 2 and sin 2 29 lying in the intervals 
4.7 x 10~ 12 eV 2 £ Am 2 £ 6.5 x 10~ 12 eV 2 , 0.71 £ sin 2 29 < 1.0, (14) 

(B) for any value of f Be from the interval (9) (for f Bo = 0.7) and 0.45 (0.42) £ f B £ 0.65 (0.66), 
and for Am 2 and sin 2 29 having values within the intervals 

3.2 x 10~ n eV 2 £ Am 2 £ 4.0 x 10~ n eV 2 , 0.59 £ sin 2 20 < 1.0. (15) 

Both solutions (A) and (B) are stable with respect to changes of f Be |. Nevertheless the 
regions of values of Am 2 and sin 2 29 of these solutions vary somewhat with f Be : eqs. (14) 
and (15) represent the largest intervals and correspond practically to f Be = 0.7. 

Let us discuss the above results. The probability that a solar electron neutrino with 
energy E will not change into &v T ) (or u s ) on its way to the Earth when v e <-> v^m [y e ^ v s ) 
oscillations take place, can be written in the form: 

P osc (E; R(t)) = 1 - I sin 2 29 [1 - cos 2t3^], (16) 
2 L v 

where L v = 47rE/Am 2 is the oscillation length in vacuum, 

R(t) = Ro [l-ecos27r^], (17) 

is the Sun-Earth distance at time t of the year (T = 365 days), Ro = 1.4966 x 10 s km and e = 
0.0167 being the mean Sun-Earth distance and the ellipticity of the Earth orbit around the 
Sun. The term with the e factor in eq. (17), as is well known [2-4,25], is a source of seasonal 



2 For f Bo = 0.7 there are also new solutions in the region 10~ 10 eV 2 < Am 2 < 10~ 9 eV 2 and 
0.62 ^ sin 2 29 ^ 0.80, representing three very narrow strips (almost lines) of allowed values 
of Am 2 and sin 2 29 in the Am 2 — sin 2 29 plane (see Fig. la). However, these solutions are 
not stable with respect to variations of f Be and disappear when f Be is slightly increased (they 
do not exist for f Be = 1.0, for example). We shall not discuss them further. 

3 The possibility of a "low" 8 B neutrino flux solution for f Be = 1.0 at Am 2 = 6.0 x 10~ 12 eV 2 
and sin 2 29 = 0.8 was suggested on the basis of qualitative arguments in ref. [29]. Our results 
show that at 95% CL. the indicated point in the relevant parameter space is marginally 
excluded by the current solar neutrino data . 
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variation effects in the case of the vacuum oscillation solution of the solar neutrino problem 
Since we are dealing in the analysis of interest with experimental results averaged over at least 
few complete years of measurements, the relevant probability is actually the probability (16) 
averaged over a period of 1 year, P osc (E; R , e). If 27rR /L v ^ 2.5n, P osc (E; R , e) practically 
coincides with the probability (16) in which the parameter e is formally set to zero, i.e., with 
P osc (E;Ro) (see Figs. 2a and 2b). This implies that for the values of Am 2 ^ 10~ 10 eV 2 of 
interest one has P osc (E; R , e) = P osc (E; Ro) for all neutrinos with energy E ^ 3 MeV, i.e., 
for the dominant fraction of the 8 B neutrino flux. If, however, 27rRo/L v >> 2.5tc, the effect 
of the averaging can be quite dramatic for the oscillation's amplitude [] and (for a given 
sin 2 26) P osc (E; R , e) can differ considerably from P osc (E; R ), as Figs. 2a and 2b illustrate. 

Consider first the solutions at Am 2 ^ 4.4 x 10 -11 eV 2 , i.e., in the region in which the 
fe = fee = 1 solution lies. The allowed regions corresponding to these solutions converge 
continuously (changing their shape and dimensions) to the allowed region in the case fe = 
fee = 1 when fe and fe e are varied continuously from the values they have for a given 
solution to 1. The new solutions (A) and (B) identified above are "disconnected" from the 
fB = fBc = 1 solution and disappear when f B and fee change continuously from 0.35 to 1. 

For Am 2 ^ 4.4 x 10 -11 eV 2 and for the energies of the pp neutrinos, E < 0.42 MeV, the 
cosine term in the expression for the probability (16) is a fastly oscillating function of E. 
Therefore the integration over the neutrino energy E in the contribution of the pp neutrinos 
to the signal in the Ga-Ge experiments suppresses the part of the contribution containing 
the cosine term and one has effectively in this case [2,3] Ppp,(E;R ) = 1 — |sin 2 26 l . This 
implies that depending on the value of sin 2 26 the pp v e flux is suppressed approximately by 



4 Detailed predictions for the seasonal variation effects in the present and future solar neu- 
trino experiments in the case of the fe — 1, fBc — 1 vacuum oscillation solution are given in 
ref. [25]. 

5 It is not difficult to convince oneself that up to corrections which do not exceed 5 x 10~ 3 the 
periods of v e <-> v^m oscillations implied by the probabilities P osc (E;R ,e) and P OSC (E;R ) 
coincide. 
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a factor (0.50 - 0.70). The values of Am 2 for which the minimal values of sin 2 29 allowed by 
the data occur (see Figs, la - lc) are determined by the condition cos27rR /L v = —1 for E 
= 0.862 MeV - the energy of the dominant 7 Be neutrino component of the solar neutrino 
flux; they correspond for a given sin 2 29 to a maximal possible suppression of the flux of the 
0.862 MeV 7 Be electron neutrinos at the Earth surface due to the oscillations u e <-> v^)- 
Actually, the suppression of the 0.862 MeV 7 Be v e flux due to the vacuum oscillations is not 
very strong and can practically be absent in the case of the solution under discussion, the 
relevant suppression factor ranging from approximately 0.30 (possible when f Be = 1.3) to 
0.98 (possible for f Be = (0.7 - 1.0) and f B = 1). 

The maximal value of f B , max f B — 3.4, for which the u e <-> v^r) oscillations provide 
a description of the solar neutrino data is determined primarily by the Kamiokande result 
(2) (as its independence on f Be indicates) and by the specific dependence of P OSC (E;R ) (see 
eq. (16)), on the solar neutrino energy E. It can be understood qualitatively by consid- 
ering the constraints the Kamiokande data imply in this particular case. For any fixed 
Am 2 ^ 10~ 10 eV 2 , the probability P OSC (E;R ) has at most one minimum in the interval of 
8 B neutrino energies 7.5 MeV ^ E < 14.4 MeV relevant to the Kamiokande experiments (see 
Figs. 2a and 2b). The suppression of the integral in the expression (6) for R(K) is maximal 
when the minimum of P osc (E;Ro) occurs at E = E min , 7.5 MeV < E min < 14.4 MeV. In 
this case P osc (E; R ) increases monotonically with the change of E in the indicated interval 
both for E < E min and E > E min . This implies that there exists a maximal possible sup- 
pression of the flux of 8 B electron neutrinos with E ^ 7.5 MeV (i.e., of the integral in the 
right hand side of eq. (6)) due to P osc (E; R ), and hence a maximal possible value of f B for 
which the vacuum v e <-> v^t) oscillations can provide an explanation of the Kamiokande 
result (2). The value one obtains numerically, max f B = 3.4, is somewhat smaller than 
the upper bound (8) and is reached for, e.g., f B = 1 at Am 2 = 8.2 x 10~ n eV 2 and 
sin 2 29 = 1.0. For these values of Am 2 and sin 2 29 one has minP osc (E min ; R ) = and 
E m i n = 9.5 MeV, and the suppression of the integral in eq. (6) corresponds effectively 
to a constant factor [P OSC (E;R ) + 0.16(1 - P OSC (E;R ))] = 0.2. The value of Am 2 for 
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which the solution for f B = 3.4 exists depends, although weakly, on the value of fee cho- 
sen within the interval (9) because P osc (E; Ro, e) for E=0.862 MeV is very sensitive to even 
small changes of Am 2 in the vicinity of 8.2 x 1CT 11 eV 2 , as Fig. 2b (after the necessary 
rescaling of the values of E on the horizontal axis by the ratio (lCT 10 /(8.2 x 1CT 11 ) = 1.2) 
illustrates (we have: 2irR /L v = 37.9018(Am 2 /l(r 10 eV 2 )(lMeV/E) and in this case, for 
instance, P osc (E = 0.862 MeV;R ,e) = P osc (E = 0.862 MeV; R ) = 0.50 (0.70) for Am 2 = 
8.2 (8.3) x 10 -11 eV 2 and sin 2 29 = 1.0). Therefore the necessary suppression of the 7 Be 
contribution to the signals in the Cl-Ar and the Ga-Ge experiments for the different val- 
ues of fee considered is achieved by small changes of the value of Am 2 around the value 
8.2 x 10~ n eV 2 . These changes do not affect the suppression of the contributions of the 8 B 
neutrinos in R(Ar) and R(K), required by the Cl-Ar and Kamiokande data. 

In a similar way one can understand the minimal values of f B , min f B = 0.51; 0.54; 0.57 
corresponding to f Be = 0.7; 1.0; 1.3, for which there exists at 95% C.L. au e ^ v^ T ) oscillation 
solution in the region Am 2 ^ 4.4 x 10~ n eV 2 . Depending on f Bc these solutions occur for 
values of Am 2 ^ (4.3-5.1) x lO" 11 eV 2 and for sin 2 29 ^ (0.57-0.75) (see Figs, la-lc). The 
minimal values of f B allowed by the data are again determined by the Kamiokande result, 
and by the specific energy dependence of P OSC (E;R ) for values of Am 2 in the vicinity of 
Am 2 = 4.9 x 10~ n eV 2 . Indeed, for Am 2 = (4.3-5.1) x 10~ n eV 2 and for the energies of 8 B 
neutrinos 7.5 MeV ^ E < 14.4 MeV detected by the Kamiokande experiments, P OS c(E; Ro) is 
a monotonically (rather steeply) increasing function of E and for sin 2 29 = 0.60 (0.75) one has 
0.5 (0.4) £ P osc (E; R ) £ 0.8 (0.7) (see Figs. 2a and 2b). In this case and, e.g., for f Be = 0.7, 
the maximal suppression due to P OSC (E;R ) of the integral in the expression for R(K), eq. 
(6), corresponds effectively to a constant factor [P OSC (E;R ) + 0.16(1 — P osc (E;Ro))] = 0.7. 
The Kamiokande data then imply (95% C.L.) f B ^ 0.52 . The minimal value of f B one 
obtains depends somewhat on the value of f Be because the requisite suppression of the 7 Be 
0.862 MeV electron neutrino flux (and of the contributions of the 0.862 MeV 7 Be neutrinos to 
R(K) and R(Ge)) is achieved now by an adjustment of the value of sin 2 29 within the interval 
0.57 - 0.75 (rather than by changing Am 2 ), which in turn leads to a non- negligible change 
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of P osc (E; R ) for 7.5 MeV £ E < 14.4 MeV. In the case of the solution with f Be = 0.7 and 
f B = 0.53, for instance, the pp and 0.862 MeV 7 Be electron neutrino fluxes are suppressed due 
to the oscillations by the factors 0.70 and 0.46, respectively, while the pep (CNO) neutrino 
flux (fluxes) is not (are mildly) suppressed. The predictions for the signals in the Cl-Ar and 
Ga-Ge detectors read in this case R(Ar) = 2.7 SNU and R(Ge) = 74 SNU. 

It should be noted that the allowed regions found at 95% C.L. for fee — (0.7 — 1.3) and 
f B = (0.8 — 1.2) in the present study lie practically all within the allowed regions one obtains 
at 95% C.L. in the reference model [7] when the estimated uncertainties in the theoretical 
predictions are included in the analysis. 

Let us add finally that the minimal values of the x 2 ~ function for the solution under 
discussion in the three cases fee = 0.7; 1.0; 1.3 respectively read 0.67; 0.51; 0.53 (for 2 d.f.) 
and take place at f B = 2.2; 2.4; 2.4 for Am 2 = (7.2; 7.4; 7.5) x KT 11 eV 2 and sin 2 29 = 1. For 
f B < 1, minx 2 = 2.9 Q and is reached for f B = 1, Am 2 = 6.1 xlCT 11 eV 2 and sin 2 29 = 0.86. In 
the case of solution (A), eq.(14), and for f Be = 0.7 (1.0), min x 2 = 4.4 (4.5) and corresponds 
to f B = 0.39, Am 2 = 5.4 (5.6) x 1CU 12 eV 2 and sin 2 26* = 1.0. The min x 2 value is somewhat 
larger for solution (B): for f Be = 1.0, for instance, one has min \ 2 = 5.3. 

3.1.2 Physical Implications of the New Low $ B Solution 

The physical implications of the v e <-> v^m oscillation solution for f Be — 1 and f B = 1 
and values of Am 2 in the interval 4.4 x 10 -11 eV 2 ^ Am 2 ^ 10~ 10 eV 2 have been extensively 
discussed in the literature (see refs. [2,25], [23] and the articles by A. Acker et al. and by V. 
Barger et al. quoted in ref. [3]). The solutions we have found in the same Am 2 region for 
f Bc 7^ 1 and f B ^ 1 lead to generically similar implications and we shall not consider them 
here. 

Of the two new solutions (A), eq. (14), and (B), eq. (15), solution (A) is more in- 
teresting phenomenologically, has a lower x 2 — value, and therefore we shall discuss only it, 



6 In contrast, min x 2 — 0.25 in the case of the MSW small mixing-angle v e — > v^m solution. 
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although rather briefly. For the values of Am 2 from the interval given in (14) one has: i) 
P OSC (E;R , e) ^ 0.97 (0.94) for E > 7.5 (5.0) MeV, ii) for neutrino energies in the vicinity 
of 0.862 MeV P osc (E;Ro,e) has its first local minimum when E decreases from values for 
which P osc (E; Ro, e) = 1 (see Figs. 2), and iii) the first local maximum of P osc (E; Ro, e) as E 
decreases below 0.862 MeV occurs in the interval 0.23 MeV ^ E ^ 0.42 MeV Q Correspond- 
ingly, if solution (A) is realized, the signals due to the 8 B neutrinos in the present (and the 
future SNO [30] and Super-Kamiokande [31]) detectors will not practically be affected by 
the vacuum v e «-> zv T ) oscillations (remember that for solution (A) fe — (0.35 — 0.44)), the 
0.862 MeV 7 Be v e flux will be suppressed by the suppression factor (0.06 - 0.26), while the 
signal due to the pp electron neutrinos in the Ga-Ge detectors will be mildly reduced by a 
suppression factor not smaller than 0.7. Thus, from the point of view of how the different 
components of the solar neutrino flux are affected, the vacuum oscillation solution (A) is 
very similar to the low fe MSW v e — > v^ir) [21] or v e — > v s [26] transition nonadiabatic 
solution. However, some of the physical implications of the two solutions differ considerably. 
In particular, i) the predicted distortion of the spectrum of the pp neutrinos is much stronger 
in the case of the vacuum oscillation solution (A) ^ (Fig. 3a) than for the corresponding 
MSW nonadiabatic solution, and ii) if solution (A) is valid, the 7 Be and pp electron neutrino 
fluxes at the Earth surface will exhibit seasonal variations which cannot take place in the 
case of the MSW solutions (see ref. [32] and the first article quoted in ref. [24]). In what 
follows we shall discuss briefly the seasonal variation effects predicted in the case of solution 
(A). 

For Am 2 < 6.5 x 10" 12 eV 2 and E > 0.233 MeV (0.217 MeV) one has: 27reR /L v < 



7 One can explain the minimal (maximal) value of fe for which solution (A) exists, the reason 
for the difference between the maximal allowed values of Am 2 (and the values themselfs) 
for a given f& e in the cases fe = 0.35 and fe = 0.38 etc., in a similar way we did it earlier, 
e.g., for the maximal value of fe allowed by the data and the corresponding values of Am 2 . 

8 It is also very different from the distortion of the pp neutrino spectrum in the case of the 
solutions with 4.4 x 10" 11 eV 2 £ Am 2 £ 10" 10 eV 2 ( see the third article quoted in ref. [2]). 



14 



0.18 (0.19) << I. Thus, in the case of solution (A) the probability P osc (E; R(t)), eq. (16), can 
be represented as a power series in the small parameter (27re(R /L v ) cos27r(t/T)). Neglecting 
all the terms smaller than 10 -3 in this series we obtain: 

P osc (E; R(t)) = P osc (E; R ) + P scas (E; R , e, t), (18) 

where the term 

1 RR tl RRt 
P scas (E;R ,e,t) = - sin 2 29 [2ne-^ sin27r-^] cos 2tt--- sin 2 29 cos2vr-^ [27re-^ cos2tt-] 2 

2 L v L v 1 4 L v L v 1 

(19) 

is responsible for the seasonal variation effects of interest. Obviously, for fixed values of 
the parameters Am 2 and sin 2 29, the difference between the values of P seas (E; R , e, t) in 
December/ January (t = 0) and June/July (t = 0.5T) is the largest. 

For 8 B neutrinos with E > 5 MeV (6.44 MeV) we have 2vreRo/L v < 8.2 x 10~ 3 (6.3 x 10~ 3 ) 
and it follows from eqs. (18) and (19) that the seasonal variation effect in the signal of the 
Super-Kamiokande (SNO) detector will be too small to be observable. The effect can be 
much larger for the signals due to the 7 Be and/or pp neutrinos in the Ga-Ge, BOREXINO 
[33] and HELLAZ [34] detectors. 

In the case of solution (A) one has for the predicted average rate of Ge production per year 
in the Ga-Ge experiments for f Be = 0.7; 1.0; 1.3: R(Ge) £ 80; 84; 88 SNU. The difference 
between the rates of Ge production in December/ January (t = 0) and June/ July (t = 0.5T), 
AR seas (Ge), due to i) the term (19) in the vacuum oscillation probability (18), and ii) the 
change of the neutrino fluxes with the change of the Sun-Earth distance due to the standard 
geometrical effect, as can be shown, satisfies: 4.0 SNU ^ AR seas (Ge) ^ 8.7 (8.1) SNU, 
the maximal value corresponding to fee = 1-3 (0.7). A convenient relative measure of the 
predicted seasonal effect is the seasonal (December/ January - June/July) asymmetry: 

A scas (Ge) - _^ L^y-No - -^-|t.o.5T], (20) 

where R(Ge;t) [Ro/R(t)] 2 is the rate of Ge production at time t of the year and R(t) is 
given by eq. (17). For solution (A) we have: 0.072 ^ A seas (Ge) ^ 0.13, the contribution due 
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purely to the geometrical factor R~ 2 (t) being 4e = 0.0668. For given Am 2 and sin 2 29 the 
change of A seas (Ge) with the change of f Be is negligibly small. 

The corresponding seasonal (December/ January - June/July) asymmetry in the signal 
due to the 7 Be (0.862 MeV) neutrinos in the BOREXINO detector is given (up to corrections 
smaller than 10~ 3 ) by 

A s a eas (BOR) = 0.0668 + °- 79 ^W^v 2 Q 2 Ro ( } 

seasV ; 0.21 + 0.79P OSC (E; R ) L v ' v ; 

where 0.0668 is the asymmetry in the absence of oscillations, and we have used the fact that 
a u ^ r)e (E)/a Uee (E) = 0.21 for E = 0.862 MeV. Note that the asymmetry A^ cas (BOR) does 
not depend on the total flux of 0.862 MeV 7 Be neutrinos. Note also that in the case of 
solution (A) for E = 0.862 MeV we have sin27rR /L v > and the second term in eq. (21) is 
always positive. As can be shown, the asymmetry Ag Cas (BOR) changes very little with the 
variation of Am 2 and sin 2 26 within the allowed regions of values for the solution (A) and 
A^ cas (BOR)- (0.11 -0.13). 

The HELLAZ experiment [34] is envisaged to detect pp neutrinos having energy 
E > 0.217 MeV and to measure their spectrum. The experiment will be based on the v — e~ 
elastic scattering reaction. Since the energy of the incident pp neutrino in each event will 
be reconstructed, one can define a seasonal asymmetry in the signal of HELLAZ, gener- 
ated by neutrinos having energy within a given interval Ei < E < E 2 (Ei > 0.217 MeV, 
E 2 < 0.42 MeV): A seas (H; E 1; E 2 ). The expression for A scas (H; E 1; E 2 ) can be obtained for- 
mally from eq. (20) by replacing R(Ge; t) and R(Ge) with the corresponding quantities - 
event rate at time t of the year, R(H; Ei, E 2 , t), and mean event rate per year, R(H; Ei, E 2 ), 
for HELLAZ. 

For solution (A) of interest and E > 0.217 MeV one has 27reRo/L v ^ 0.19, and the sea- 
sonal asymmetry in the total (neutrino energy integrated) signal of the HELLAZ detector, 
Ag eas (H), as numerical calculations show, satisfies 0.018 ^ Ag Cas (H) ^ 0.067, the smallest 
(the largest) value being reached for Am 2 = 5 x 10~ 12 eV 2 (6 x 10" 12 eV 2 ). Note that for 
Am 2 = 5 x 10~ 12 eV 2 the asymmetry due to the v e <-> oscillations compensates to a large 
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extent the geometrical one, rendering the net asymmetry hardly observable. The absence 
of any seasonal variations in the signal of SNO, or Super-Kamiokande, or BOREXINO, or 
HELLAZ detector (constant in time event rate) is known [2,25] to be one of the distinctive 
signatures of the vacuum oscillation solutions of the solar neutrino problem. 

The comparatively small values of the asymmetry Ag Cas (H) are a consequence of two 
circumstances. First, the contribution of the v^m neutrinos (present in the pp (u e ) flux as 
a result of the oscillations) reduces the asymmetry of interest. For a given energy E the 
expression for the latter contains the factor (1 — er„ ,. e (E)/ ov ee (E)) which changes from 0.58 
to 0.72 when E increases from 0.217 MeV to 0.42 MeV []. More importantly, for values of 
Am 2 from the interval (14), sin27rR /L v changes sign passing through zero in the interval 
0.28 MeV < E < 0.39 MeV when E varies from 0.217 MeV to 0.42 MeV 0. As a result 
the two contributions to the asymmetry in the energy integrated event rate, generated by 
the oscillations of pp neutrinos having energy in the two intervals 0.217 MeV < E < Eo and 
E < E < 0.42 MeV, E > 0.28 MeV being the energy at which sin2vrR /L v = Q, have 
opposite signs and compensate partially or completely each other. A complete cancelation 
between the indicated two contributions takes place, for instance, for Am 2 = 6 x 10 -12 eV 2 , 
for which E = 0.36 MeV. 

It should be evident from the above discussion that for solution (A) the asymmetry 
A: eas (H; 0.217 MeV, E ) = A* eas (H; E < E Q ) or A s a eafi (H; E , 0.42 MeV) = A^ as (H;E>E ) 
can be larger than Ag Cas (H). Indeed, it can be easily shown that either |Ag eas (H; E < E )| > 
A^ cas (H), or |A- eas (H;E>E )| > A s a eas (H). One can have also |A* cas (H; E < E ) - 



9 The same term gives rise to the factor 0.79 in the expression for Ag eas (BOR), eq. (21). 

10 For 5.4 x 10~ 12 eV 2 ^ Am 2 ^ 6.5 x 10~ 12 eV 2 , sin27rR /L v changes sign two times 
in the interval (0.217 MeV - 0.42 MeV), passing through a second zero located at 
0.217 MeV < E < 0.26 MeV . However, the effect of the presence of this second zero 
of sin27rR /L v on Ag eas (H) is less important than the effect of the first zero located at 
E > 0.28 MeV. 

11 Obviously, the value of E depends on the value chosen of Am 2 . 
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A a cas (H; E > E )| > A s a eas (H). Note that the difference [A s a cas (H; E < E ) — A s a cas (H; E > E )] 
is free from the geometrical term and can be nonzero only if the pp neutrinos take part 
in vacuum oscillations. Let us illustrate the above remarks with two examples. For 
Am 2 = 6 x 1(T 12 eV 2 (E = 0.36 MeV) and sin 2 29 = 1 one has: A a cas (H) = 0.0668 and 
A a cas (H; E < 0.36 MeV) 0.11. If Am 2 = 5 x 10~ 12 eV 2 , then E = 0.30 MeV and in this 
case A a cas (H) = 0.018, A s a cas (H; E < 0.30 MeV) 0.10, A a cas (H; E > 0.30 MeV) £ lO" 3 , and 
[A a eas (H; E < 0.30 MeV) - A s a cas (H; E > 0.30 MeV)] = 0.10. 

In the case of the solution with 4.4 x 10~ n eV 2 ^ Am 2 ^ 10~ 10 eV 2 the seasonal asym- 
metry in the signals of the Ga-Ge, Super-Kamiokande (SNO) and BOREXINO detectors 
due purely to the v e <-> u^ T ) oscillations can be as large as 30%, 14% and 80%, respectively, 
and is predicted to be negligible for the energy integrated signal of the HELLAZ detector 
[2,25]. 

3.2 Oscillations into Sterile Neutrino v e <-> v s 

The solar neutrino oscillations into sterile neutrino, u e <-> u s , were excluded in the case 
f B = 1 and f Be = 1 at 99% C.L. as a possible solution of the solar neutrino problem by 
the (mean event rate) solar neutrino data which were available by March 1994 [4]. Since 
then updated results have been published by all operating solar neutrino experiments. The 
current status of the hypothesis of solar u e <-> v s oscillations, including the cases f B ^ 1 and 
f Be 7^ 1, 0.7 < f Be < 1.3, is summarized graphically in Figs. 4a - 4c. At 95% C.L. and for 
0.7 < f Bc < 1.3 (f Bc = 0.7) this possibility is not excluded by the current solar neutrino data 
only for 

0.35 £ f B £ 0.43 (0.44) (22) 
and values of Am 2 and sin 2 29 in the intervals 

4.8 x 1(T 12 eV 2 £ Am 2 £ 6.2 x 1(T 12 eV 2 , 0.74 £ sin 2 29 < 1.0. (23) 

This solution is stable with respect to variations of f Bc within the interval (9). Obviously, 
it is a u e <-> v s oscillation analog of the v e <-> v^) oscillation solution (A) (compare eqs. 
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(22), (23) with eq. (14)). At 95% C.L. and for f Be = 0.7 there exist also two other allowed 
regions at larger values of Am 2 for f B = 0.42 and f B = 0.50 (see Fig. 4a), but they are very 
small and disappear when f Be > 0.7. If one increases the required C.L. of the description of 
the data to 98%, solution exists for 0.32 ^ f B ^ 1.8 if f Bc = (0.7 — 1.3). The corresponding 
new allowed regions are scattered over the area 4.0 x 1CT 12 eV 2 ^ Am 2 ^ 2.0 x 1CT 10 eV 2 , 
0.50 ^ sin 2 29 < 1.0. These regions diminish in size considerably or completely disappear 
as f Be changes from 0.7 to 1.3: for f Be = 1.3 most of the remaining new regions are in the 
form of narrow strips (see Fig. 4c). 

Let us consider briefly the physical implications of the solution (22) - (23) for the future 
solar neutrino experiments. The deformation of the pp neutrino spectrum (Fig. 3b) is quite 
strong and differs somewhat from the deformation in the case of the v e <-> z/ M(r ) oscillation 
solution (A) (compare Figs. 3a and 3b). In what regards the seasonal variation effects in 
the signals of the Super-Kamiokande, SNO and the Ga-Gc experiments, they coincide with 
those for the v e <-> oscillation solution (A), considered in Section 3.1.2. The predicted 
seasonal variation effect in the signal due to the 7 Be neutrinos in the BOREXINO detector, 
however, is larger than in the corresponding case of v e <-> v^r) oscillations. The seasonal 
(December /January - June/ July) asymmetry for the BOREXINO detector is given now by 
the expression 

ALs(BOR) = 0.0668 + sin 2 29 sm2n^, (24) 

and can be as large as 42%: one has 0.18 ^ Ag eas (BOR) ^ 0.42. We find also that 
A^ cas (BOR) k, 1.5A s a cas (BOR). The difference between the values of A^ cas (BOR) and 
Ag Cas (BOR), in particular, can be used to distinguish between the u e <-> z/ M(r ) oscillation 
solution (A) and its v e <-> v s oscillation analog in a solar model independent way. 

The seasonal asymmetries Ag eas (H;E< E ) and Ag eas (H;E> E ) in the signal of the 
HELLAZ detector tend also to be larger than the corresponding asymmetries in the case of 
the v e <-> z/ m(t) oscillation solution (A). For Am 2 = 5.0 x 1(T 12 eV 2 (E = 0.30 MeV) and 
sin 2 29 = 1.0, for instance, we have: A^ cas (H) -0.9 x 10^ 2 , A^ cas (H; E < 0.30 MeV) = 0.12, 
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A^ cas (H; E > 0.30 MeV) - -0.04, and [A* cas (H; E < 0.30 MeV) - A* eas (H; E > 0.30 MeV)] - 
0.16. In the case of Am 2 = 6.0 x 10~ 12 eV 2 (E = 0.36 MeV) one obtains: A* cas (H) = 0.0668, 
and A^ cas (H; E < 0.36 MeV) = 0.15. 

4. ENERGY INDEPENDENT SUPPRESSION OF THE SOLAR 
NEUTRINO FLUX 

The possibility of universal (energy independent) suppression of the pp, 7 Be, pep, 8 B and 
CNO neutrino fluxes can be realized if solar neutrinos take part in v e <-> or u e <-> u s 
oscillations characterized by Am 2 >> 10~ 4 eV 2 . The solar matter effects for Am 2 >> 
10" 4 eV 2 are negligible and neutrinos propagate in the Sun as in vacuum. The averaging 
over the region of neutrino production, etc. in the indicated case renders the oscillating term 
in the expression for the oscillation probability, eq. (16), negligible and one effectively has 
Pose = 1 — 1/2 sin 2 29 for all components of the solar neutrino flux. The Voloshin, Vysotsky, 
Okun [35] solar u e spin precession scenario also leads to the indicated type of reduction of 
the solar v e flux. 

In general one has to consider two possibilities: transitions (or oscillations) into active 
neutrino, v e — > v^t) or u e — > P M ( T ), and into sterile neutrino u e — > v s . From the point of view 
of the analysis of the solar neutrino data currently available, there is no difference between 
the cases of u e — > v^r) and v e — > transitions (or oscillations). This follows from the 
fact that for E ^ 7.5 MeV the cross-sections <7^ (r)e (E) and cr p ^ (T)e (E) practically coincide. 

We have investigated the possibility that the solar neutrino deficit is due to a suppression 
of the different components of the solar neutrino flux by one and the same energy independent 
factor R resulting from u e — > [y e <-> z^( T )) or v e — > z^( T ), or from v e — > v s {y e <-> u s ) 

transitions (oscillations). There are two parameters in the corresponding x 2 -analysis: R 
and f B . They were varied within the intervals: (0.0 - 1.0) and (0.0 - 5.0), respectively. The 
parameter f Bc was assumed to have a fixed value within the interval (9). 

Our analysis showed that for fe e = 0.7; 1.0; 1.3 a neutrino energy independent suppression 
of the solar neutrino flux resulting from u e — > v^r) { v e ^ u ^(t)) or v e — > transitions 
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(oscillations) is excluded by the current solar neutrino data at 97%; 98%; 98% C.L. The 
regions in the R - f B plane allowed in this case at 99% C.L. (x 2 < 9.21) are shown in 
Figs. 5a - 5c. Finally, for the indicated values of f Be the solar neutrino data rule out the 
hypothesis of constant suppression of the solar neutrino flux due to u e — > u s [y e <-> u s ) 
transitions (oscillations) at 99.0%; 99.5%; 99.7% C.L. 

5. CONCLUSIONS 

We have shown that the u e <-> v^t) vacuum oscillation solution of the solar neutrino 
problem is stable with respect to changes in the predictions for the fluxes of 8 B and 7 Be 
neutrinos. For low values of $ B (fB — 0.35 — 0.43) new u e <-> u^ T ) and v e <-> v s oscillation 
solutions exist. We have discussed the physical implications of these new solutions for the fu- 
ture solar neutrino experiments. A second new v e <-> v^{ T ) oscillation solution has been found 
for values of f B which lie within the interval f B = 0.45 — 0.65. The current solar neutrino 
data exclude at 99 % C.L. the possibility of universal (energy independent) suppression of 
the different components of the solar neutrino flux, caused by u e <-> v s oscillations or u e — > v s 
transitions. A similar suppression resulting from solar u e oscillations or transitions into an 
active neutrino [y e <-> z^( T ), v e — > z^( T )) is strongly disfavoured by the data: depending on 
the value of f Be it is excluded at 97%-98% C.L. 
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Figure Captions 

Figs, la - lc. Regions of values of the parameters Am 2 and sin 2 29 for which the solar 
neutrino data can be described at 95% C.L. in terms of u e <-> v^t) oscillations of the solar u e 
for values of f Be = 0.7 (a), 1.0 (b), 1.3 (c), and for values of f B from the interval (0.35 - 2.5). 

Figs. 2a - 2b. The vacuum oscillation probability for the mean distance between the 
Sun and the Earth, P OSC (E;R ) (a), and the probability (16) averaged over a period of 1 
year, P osc (E;R ,e) (b), as function of the neutrino energy E for Am 2 = 10~ 10 eV 2 and 
sin 2 29 = 0.8. 

Figs. 3a — 3b. The deformation of the normalized to one spectrum of pp neutrinos in 
the cases of a) v e <-> v^t) an d b) u e <-> v s oscillation solutions (14) and (23), respectively. 
The dotted, dashed, long-dashed, dash-dotted and long-dash-dotted lines correspond to 
Am 2 = (5.2; 5.4; 5.6; 5.8; 6.0) x 10~ 12 eV 2 and sin 2 29 = 1. 

Figs. 4a - 4c. The regions of values of the parameters Am 2 and sin 2 29 for which the solar 
neutrino data can be described at 95% C.L. (dashed lines) and 98% C.L. (solid lines) in 
terms of u e <-> v s oscillations of the solar v e for f Be = 0.7 (a); 1.0 (b); 1.3 (c), and for values 
of f B from the interval (0.35 - 1.5). 

Figs. 5a - 5c. The regions of values of the parameters R and f B allowed at 99% C.L. (% 2 < 
9.21) by the solar neutrino data in the case of universal (energy independent) suppression 
of the different components of the solar neutrino flux by one and the same factor R, caused 
by vacuum oscillations or transitions of the solar neutrinos into an active neutrino (^ M ( T ) or 
*V( T )) in the three cases f Be = 0.7 (a); 1.0 (b); 1.3 (c). 
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